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ABSTRACT 

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an on-going Hubble Space Tele- 
scope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 
disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment 
of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the 
largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial 
resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously 
inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, repre- 
senting more than a factor of four increase over previous catalogs within the current survey area (390 
arcmin 2 ). This work presents results derived from the first ~25% of the survey data; we estimate that 
the final sample will include ~2500 clusters. For the Year 1 objects, we present a catalog with posi- 
tions, radii, and six-band integrated photometry. Along with a general characterization of the cluster 
luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and 
highlight a number of individually interesting clusters found in the Year 1 search. 
Subject headings: catalogs — galaxies: individual (M31) — galaxies: star clusters: general 



1. INTRODUCTION 

Large, high-quality samples of stellar clusters provide 
key data for studies of a wide variety of astrophysical 
topics, including cluster evolution, stellar evolution, star 
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formation, and galaxy evolution. However, current clus- 
ter samples suffer from serious limitations. For example, 
Milky Way clusters suffer from severe dust attenuation 
within the Galactic plane, resulting in disk cluster sam- 
ples that are complete only within a small region (ra- 



Piskunov et al.|[2QQ8). This limits the variety of objects 



dius of ^1 kpc) arou nd the Sun (e.g., |Dias et aT][2QQ2 

obi 



and galactic environments explored, as exemplified by 
the dearth of massive, inter mediate-age clusters k nown 



within the Milky Way (e.g., Davies et al. 2011a). Al- 



though infrared surveys of the Galactic plane are improv- 
ing this situation (e.g., putra et al.||2QQ3| |Mercer et al.| 



2005] |Borissova et al. 



201 ip , current samples of Galactic 
clusters do not probe the full stellar cluster parameter 
space, limiting our ability to study mass, age, and envi- 
ronmental dependencies of evolutionary processes. 

One solution to the incompleteness of Galactic samples 
has been to extend cluster studies to other galaxies. Ex- 
tragalactic cluster samples have grown immensely over 
the past decade, thanks in part to the power of Hubble 
Space Telescope (HST) imaging. These analyses probe 
a variety of galactic environ ments, ranging from st ar- 
bursting galaxy merg ers (e.g.. I Whit more et al.||1999[ ) to 
quiescent spirals (e.g., |Larsen||2002| ), producing provoca- 
tive results concerning the environmental dependence of 



cluster formation and evo lution (e.g., Larsen & Richtler 
20001 IGoddard et al.||2010|). However, even with HST 



clusters in most distant galaxies appear as marginally 
resolved single objects. This increases the difficulty and 
uncertainty associated with basic steps in cluster anal- 
ysis, such as object identification, photometry, and the 
derivation of ages and masses. As a result, the interpreta- 
tion of underlying cluster evolutionary processes has con- 
siderable associated uncertainty, even leading to cases of 
conflicting interpretations derived from the same dataset 
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(e.g., in M83; IChandar et al.| 2010b| |Bastian et aT]|2Q12 
Fouesneau et al.|2012|) . 

Closer to the Milky Way, large cluster catalogs exist for 
the Large and Small Magellanic Clouds (LMC fc SMC; 
Bica et al.||2QQ8| [1999) |Bica fc Dutra||2QQ01 |Hunter et al. 



2003 ) . In the Clouds, the ability to resolve clusters into 
individual stars has resulted in a number of important re- 
sults in ste l lar evolution (e.g., jChiosi et al.||1989 Frogel 



et al . 1990; Gira rdi et al.||2009|), cluster evol ution (e.g. 
Gieles fc Bastian 2008; Chandar et al. 2010a), and clus- 



ter formation (e.g., |lVlackey et al.||2008| ). However, there 
are limitations associated even with this excellent sample 
of objects. On a practical level, while these cluster sam- 
ples do not suffer from the same incompleteness issues as 
those from the Milky Way, their piecemeal assembly by 
many different groups has imprinted a complex and little- 
understood selection function. On a more fundamental 
level, there are additional limits to the applicability of 
stellar and cluster evolution results derived from these 
interacting, relatively low-mass galaxies. The LMC and 
SMC do not probe the range in galactic environments 
that are characteristic of a majority of the baryonic uni- 
verse. More than 75% of all stars in the Universe today 
have metallicities within a factor of two of the solar value 



(Gallazzi et al. 2008), higher than those probed by the 
Clouds. JNo cluster sample comparable to those in the 
Clouds is currently available in a large spiral galaxy. 

The neighboring galaxy M31 (the Andromeda galaxy) 
is a prime target for stellar cluster studies. The galaxy's 
proximity allows resolution of individual bright stars in 
clusters and the robust detection of clusters down to faint 
(< 10 4 L ) luminosities. M31 also provides access to a 
range of widely-varying environments across the extent 
of the star-forming disk. Andromeda's role as a valuable 
laboratory was realized long ago, and decades of work 
have gone into exploring its cluster population. We defer 
the detailed review of exis ting clus ter catalogs until Sec- 
tion [5j but beginning with |Hubble| (|1932 ) , previous work 
has mainly focused on M3rs globula r cluster population 
through ground-bas ed imaging (e.g., |Galleti et al7| 2004; 
Peacock et al.||2010[ ). The proximity ol M31 also enables 
high- resolution spectroscopic follow-up of b right clusters 
(e.g., |Colucci et al.|[2009| jStrader et al.||2011[ ), as well as 
low-resolution spectroscopy ol intrinsically fain t clusters 
that are in accessible in distant galaxies (e.g., |Cal 
et al.| |2009|). Numerous studies have utilized H5T 
cellen 



Caldwell 



s ex- 



2009k. Numerous studies have utilized 
spatial resolution to study massive clu sters (e.g., 



Barmby fc Huchra||2001 
individual resolved stars 



Perina et al. 2009b) and their 



(e.g., 



Rich et al. 2005; |Mackey 



et al.||2006||Perina et al.|2009a|^UTT| , as welTas to iden- 
tity and study less massive disk clusters (e.g., the Hodge- 
Krienke Catalogs ; |Krienke fc Hodge|[2007l [20081 |Hodge| 
eT!I1[2009l [20101 hereafter the HKC). While extremely 
valuable, these previous space-based observing programs 
have been limited to focused studies of a small number 
of targeted regions, as opposed to a wide-ranging survey 
to obtain a broad sampling of the galaxy and an overall 
assessment of the M31 cluster population. 
The Panchromatic Hubble Andromeda Treasury 



(PHAT; palcanton et al. 
ize the study of stellar c. 



2012[ ) is poised to revolution- 
This on-going 



lusters in M31. 
HST multi-cycle program will image one-third of the M31 
disk at high spatial resolution, with wavelength coverage 
from the ultraviolet through the near-infrared. In terms 



of cluster studies, this survey provides a number of dis- 
tinct advantages over existing work. High spatial resolu- 
tion imaging allows us to resolve clusters into individual 
stars, permitting detailed characterization of their stel- 
lar populations through analysis of their color-magnitude 
diagrams (CMDs). Figure [^demonstrates the data qual- 
ity provided by PHAT for a previously unidentified clus- 
ter, showing the considerable gain over existing ground- 
based surveys. The wide wavelength coverage enables 
accurate age-dating, as well as the ability to probe a 
broad range of stellar effective temperatures, from mas- 
sive main sequence stars to evolved supergiants and AGB 
stars. As a result, these objects provide a wealth of valu- 
able knowledge in terms of calibrating and refining stellar 
evolution models at high metallicity. Finally, these im- 
ages allow the detection of faint clusters, reliably probing 
more than ~2 magnitudes further down the luminosity 
function than was previously possible using ground-based 
datasets. The PHAT survey will increase the HST spa- 
tial coverage of the M31 searched for clusters by a factor 
of ~3 (1800 arcmin 2 in total compared to 650 arcmin 2 
surveyed in the HKC; 390 arcmin 2 in the current Year 1 
dataset). However, this metric underestimates the scien- 
tific gain provided by high-quality, uniform PHAT data 
products, as opposed to the heterogeneous archival data 
used in the HKC work. Fundamentally, the PHAT sur- 
vey represents the shift from a discrete, targeted mode of 
cluster study to a broad survey mode, allowing for com- 
prehensive analysis of cluster evolutionary processes its 
environmental dependencies. 

The stellar clusters identified as part of the PHAT sur- 
vey will constitute the most comprehensive sample of 
clusters available for any large spiral galaxy. The large 
range of galactocentric radius (0-20 kpc) included in the 
survey spans a wide range of star formation intensities 
and gas densities. The diversity of galactic environments 
will be important for testing models of cluster formation 
and evolution. In addition, the simultaneous accessibility 
of objects over a >3 order of magnitude range in cluster 
luminosity provides a top-to-bottom view of the cluster 
population, given that we sample a continuous range of 
objects that extend from those equivalent to Galactic 
open clusters up to massive globular clusters. 

This paper is the first in a series utilizing the PHAT 
dataset for studies of stellar clusters. Here, we present 
the first installment of a HST-based cluster catalog, 
which will serve as the basis for extensive study of An- 
dromeda's cluster population. Catalog updates and im- 
provements will be published over the course of this 
four year observing program. In this edition, we pub- 
lish positions, sizes, and integrated photometry for the 
Year 1 cluster sample. Age and mass determinations de- 
rived from the integrated photometry will be presented 
in Fouesneau et al. (2012, in prep.). Additional stud- 
ies, including analysis of structural parameters, resolved 
star content, and integrated spectroscopy of the cluster 
sample will follow in subsequent work. 

We summarize the PHAT observations in Section [2j 
while in Section [3] we describe our cluster identification 
procedures, present results from completeness testing, 
and introduce the Year 1 cluster catalog. Next, we de- 
scribe and test our photometry methodology in Section 
[4j followed by a comparison between the PHAT cluster 
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Fig. 1. — PHAT survey data quality example, showing newly-identified cluster PC1152. Six grayscale single-band images, as well as a 
color optical (F475W+F814W) mosaic , show the superior image quality provided by HST when compared to ground-based observations 
(from the Local Group Galaxy Survey; [Massey et al.|2006| >. We also present a color-magnitude diagram of resolved photometry for objects 
that lie within the cluster cutout image. The cluster s mam sequence forms the vertical sequence at F475W-F814W~0.5, along with three 
evolved giants at F814W<20. Other stars shown with F814W>20 and F475W-F814W>1.0 are likely coincident stars belonging to the 
background field. 

ing each epoch, imaging is obtained by the cameras in 
two side-by-side, half-brick (3x3) arrays. Between the 
epochs, the orientation of the cameras change by 180 de- 
grees due to the annual roll angle variation of HST. As a 
result, the half-brick mosaic obtained by the ACS cam- 
era during the first epoch is now imaged by the WFC3 
camera in the second epoch, and vice-versa, completing 
six- filter imaging across the brick. In all, we obtain ^130 
minutes of exposure time at each of the brick's 18 field 
centers. 

The Year 1 imaging used in this work includes four 
full bricks (designated B01, B09, B15, and B21) and 
the western halves of two additional bricks (B17W and 
B23W). These data sample locations along the major 
axis of M31 from the center out to a projected radius of 
~20 kpc. The full PHAT survey footprint, along with 
the locations of the Year 1 bricks are presented in Figure 



catalog and existing catalogs in Section [5] We present 
a basic characterization of the cluster catalog contents 
in Section [6j followed by discussions of luminosity func- 
tions, the cluster size distribution, and objects of interest 
in Section [7j We conclude with a summary and descrip- 
tion of future work in Section [5J Throughout this work, 
we assume a distance mod ulus for M31 of 24.47 (785 kpc; 
McConnachie et aL]|2QQ5[ ), for which 1" corresponds to a 
physical size ol 3.81 pc. 

2. OBSERVATIONAL DATA 

This paper includes clusters identified in Year 1 PHAT 
imaging data, taken before June 2011. A full de scription 
of the PHAT ob servational design is available in|Dalcan- 
ton et al. (2012), but we briefly summarize relevant de- 
tails below. PHAT observations are grouped into 23 area 
units known as "bricks", each made up of 18 mosaiced 
HST fields of view in a rectangular 6x3 arrangement 
that covers a ~12 / x6.5 / region of sky. Data are obtained 
simultaneously at different field centers with the ACS 
(using the F475W and F814W filters) and WFC3 (using 
the F275W, F336W, F110W, and F160W filters) instru- 
ments in two epochs, separated by ~6 months. Dur- 



m 

Combining imaging data across multiple pointings 
from three separate cameras (ACS/WFC, WFC3/UVIS, 
and WFC3/IR) requires astrometry with higher preci- 
sion than that obtained using the standard telescope 
telemetry and data processing pipeline. Images of neigh- 
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Fig. 2. — Footprint of the PHAT survey region (magenta) displayed on a GALEX NUV image of the northeast half of M31. Green 
rectangles represent the "bricks" that make up the Year 1 imaging data. Blue circles show the spatial distribution of clusters identified in 
the Year 1 cluster search. 

As shown in Fig. [T] clusters appear in PHAT imaging as 
composite objects composed of a centrally concentrated 
overdensity of individual resolved stars and a broad un- 
resolved light component. The ability to resolve these 
objects into individual stars allows for the clean separa- 
tion between genuine stellar clusters and contaminants 
such as background galaxies or single stars. 

Although PHAT imaging facilitates robust identifica- 
tions of stellar overdensities, determining whether or not 
these clusterings are gravitationally bound is challeng- 
ing. Quantitative assessment of an object's bounded- 
ness requires age, mass, and spatial p rofile information 



boring fields have sufficient overlap to enable us to de- 
rive an astrometric solution across a full (or half) brick. 
These astrometric solutions are obtained separately for 
each camera using photometric catalogs deriv ed using 
DOLPHOTpJ a modified version of HSTPhot ( |Dolphin 
2000) that has been updated to include specialized ACS 
and WFC3 modules. Affine distortion corrections (in ad- 
dition to those already known for each camera from the 
IDCTABs) are required to obtain consistent brick- wide 
astrometric solutions. An additional correction brings 
the brick-wide astrometric solutions onto a global astro- 
metric frame, def ined using CFHT obs ervations tied to 
2MASS catalogs flSkrutskie et al.||2QQ6| >. The global as- 
trometric alignment agrees with that of the 2MASS ref- 
erence system within an absolute level of ~60 mas. 

Once aligned, we use the mu ltidrizzle task within 
PyRAF ( |Koekemoer et al.|2002[ ) for cosmic ray rejection 
and i mage creation along with 1 acosmic (van Dokkum 
2001[ ) for supplemental cosmic ray flagging. We note that 
cosmic ray correction is particularly difficult in the case 
of the WFC3/UVIS data, due to a large number of cos- 
mic ray artifacts (particularly in F275W images) and the 
availability of only two frames for artifact detection in 
regions of non-overlapping field coverage. Pixel scales 
of the resulting images are (0.04, 0.05, and 0.065 arc- 
sec/pixel) for the (WFC3/UVIS, ACS/WFC, WFC3/IR) 
cameras, where the WFC3/IR images are up-sampled 
from their native plate scale (0.128 arcsec/pixel) to take 
advantage of the higher effective resolution afforded by 
the survey's sub-pixel dither strategy. 

3. CLUSTER IDENTIFICATION 

The goal of the PHAT cluster survey is to identify and 
analyze a sample of gravitationally bound stars clusters 
in M31. In this paper, we undertake the first step to- 
ward achieving this goal: the visual inspection of Year 
1 PHAT imaging to identify candidate bound clusters. 



(e.g., Gieles & Portegies Zwart 2011). Further, deter 



18 http : / / purcell . as . arizona . edu/ dolphot 



mining the boundedness ol young objects (<10 Myr) is 
made even more difficult, because dynamical evolution 
has had little time to evolve stellar structures from an 
initial hierarchical, scale-f ree spatial distributi on (for fur- 
ther discussion, see e.g., Bast ian et al.||2Qll| and refer- 
ences therein), blurring the distinction between bound 
and unbound stellar groupings. In future work, we will 
utilize age, mass, and structural characteristics to assess 
boundedness for each object (e .g., using the n statistic; 
|Gieles fc Portegies Zwart |[2011| ), but that analysis is be- 
yond the scope oi the initial classification work presented 
here. 

Given the current limitations in assessing an object's 
boundedness, we adopt a liberal approach for cluster can- 
didate identification. We prioritize sample completeness 
over purity and therefore include all cluster-like stellar 
overdensities as part of the object catalog. We acknowl- 
edge that the cluster catalog presented here will likely 
include both bound and unbound groupings of stars, 
particularly among the youngest objects. As a result, 
throughout this work all objects are formally considered 
cluster candidates, although for brevity we will refer to 
them simply as clusters. 

With the intent of the cluster search well defined, 
we proceed to a description of our search methodology 
and a presentation of the cluster search results in Sec- 



5 



tion 13.11 We conclude our discussion of cluster identifi- 
cation by characterizing the completeness characteristics 
of the sample in Section [3^2) 



3.1. By- Eye Search 

We undertake a systematic by-eye search of the Year 1 
PHAT images follow ing a precedent set by previous M31 
cluster studi es (e.g., Barmby fc Huchra||2001 ; Krienke & 
Hodge 2007). Drawing on this rich history ol visual clus- 
ter identification, we make a concerted effort to improve 
upon previous work through the use of uniform analysis 
techniques, redundancy, cross-validation, and improved 
characterization of the resulting selection function. 

When using any method of cluster identification, man- 
ual or automated, it is important to understand the bi- 
ases and limitations inherent to the technique. We use 
artificial cluster tests in order to assess the complete- 
ness characteristics of our se arch methodology, which we 
discuss in detail in Section 3.2 We are also develop- 
ing automated methods of cluster identification to fur- 
ther reduce subjectivity of the identifications and enable 
more robust completeness testing in future work with this 
dataset (Olsen et al. 2012, in prep.). The cluster sam- 
ple presented here will act as an important comparison 
sample to help refine these automated search techniques 
for use with the PHAT dataset. 

Our by-eye search consists of two stages: an initial 
search of the available imaging for all viable cluster can- 
didates, followed by a re-evaluation of each preliminary 
candidate in a systematic manner. The completeness and 
accuracy of this process are enhanced by the redundancy 
of eight experienced astronomers conducting each stage 
of the search. 

To perform the initial image search, the survey area 
is subdivided using the footprint of the WFC3-IR cam- 
era (2.3 / x2.1 / ), providing 18 contiguous, minimally- 
overlapping search fields within each brick. Each field 
is searched by three or four team members in a "blind" 
manner, meaning that these individuals are not provided 
with the locations of clusters previously identified by 
other PHAT team members or previous surveys. To 
identify cluster candidates, searchers use a suite of im- 
ages that include a two-band optical color image, all six 
available single-band images spanning from the UV to 
the NIR, and two star-subtracted optical images used to 
identify diffuse emission that makes up a cluster's unre- 
solved light. These images are accessed using a custom 
image viewer that allows users to switch between the 
available spatially-aligned, full resolution images, as well 
as the ability to alter brightness and contrast levels for 
optimal visualization. 

Once candidate objects are selected in the initial visual 
search, we cross-match and combine the identifications of 
all team members. Based on this preliminary catalog, we 
perform initial photometry and create cutout images for 
each cluster candidate. The preliminary catalog is re- 
viewed independently by all eight team members on an 
interactive web site, and each individual assigns scores 
for each candidate, providing an assessment of the like- 
lihood that an object is a cluster. Scores (Sb y - eye ) are 
based on a scale of l-to-3, where Sb y - e ye = 1 represents 
a definite cluster, Sb y - eye = 2 signifies a likely cluster, 
and S b y - eye = 3 represents an unlikely or non-cluster ob- 
ject. We present examples of the cluster scoring system 



in Fig. [3j 

We average the scores from all team members for each 
object and use these average rankings to divide the can- 
didates into three subsamples: clusters, possible clus- 
ters, and unlikely objects. We choose thresholds of 



< 2.0 for clusters, 2.0 < S h . 



< 2.5 for possi- 



ble clusters, and Sb y - e ye > 2.5 for unlikely objects. We 
discard unlikely objects from the catalog, while retain- 
ing clusters and possible clusters in two separate cata- 
logs. The average scores for each cluster are provided 
as part of the cluster catalogs as an assessment of candi- 
date quality. To determine the reliability of these average 
scores, we built a ranking experiment into our classifica- 
tion procedure. During the course of the ranking work, 
24 objects appeared twice within the preliminary catalog. 
As a result, these clusters were each ranked two separate 
times by each team member. When the resulting average 
scores of the duplicate entries are compared, we find the 
standard deviation of the l-to-3 ranking differences to be 
0.27, showing good consistency and repeatability for the 
scores provided by our search team. 

The Year 1 cluster search yielded a catalog of 601 high- 
scoring clusters. Table [T] presents positions for each ob- 
ject, as well as other descriptive information (radii, pho- 
tometric measurements) that will be described in Sec- 
tion |4[ Cutout images for each object are presented in 
Fig.[4[ In addition, the spatial distribution of the clusters 
are shown in Fig. [2] Tabulated information and image 
cutouts for 237 possible clusters are presented in Table [2] 
and Fig. [5j respectively. Information about the possible 
clusters is provided for completeness, but due to the un- 
certain nature of their classifications, we exclude these 
objects from further analysis. 

During the cluster search, we also identified 370 puta- 
tive background galaxies. We did not explicitly search for 
these objects, therefore this catalog does not constitute a 
complete sample of objects. However, the potential use- 
fulness of these identifications (e.g., as astrometric refer- 
ences, multi- wavelength source catalog cross-correlation) 
warrants its release. We present this catalog of objects 
in Appendix [A] 

3.2. Catalog Completeness 

To characterize the completeness of our cluster sample, 
we conduct artificial cluster tests that mimic the selec- 
tion procedure of the clusters. We use artificial clusters 
that span the range of ages and masses we expect to find 
in the cluster sample, while the size distribution is chosen 
to sample the minimum, average, and maximum sizes of 
the true sample (see Section 7.2). Ages and masses are 
chosen from a logarithmic grid of values, while sizes are 
drawn from the set of three characteristic values. Specif- 
ically, we select ages ranging from 4 million to 10 billion 
years, masses ranging from 10 2 to 10 5 M , and profiles 
that have effective radii (Reff, equivalently, half-light 
radii) of 1, 3 or 7 pc (0.26, 0.79, or 1.84 arcsec). 

We crea te artificial cluster s by populating a Padova 
isoch rone ([Gi rardi et al.|2010| ) of the appropriate age us- 
ing a Kroupa (20(JT|) stellar initial mass function. Next, 
the st ars are spatially distributed according to a King 
(1962) profile, assuming no mass segregation. Finally, 
the size of the cluster and magnitudes of the individual 
stars are scaled appropriately to account for the distance 
of M31, assuming zero Galactic foreground or other in- 





Fig. 3. — Cluster cutout images showing objects across the range of possible Sb y - eye . The top row shows color optical (F475W+F814W) 
images, and the bottom row shows grayscale F475W images. Labels present the object name along with the average Sby-eye, which runs 
on a scale of 1 (definite cluster) to 3 (unlikely cluster). 
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Fig. 4. — Cutout images for clusters in Table [T] with St y - e ye < 
2.0. Images are grayscale F475W images, scafea to three times 
the cluster radius, and aligned such that North is up and East is 
left. Along with the PHAT cluster identifier, the average Sb y _ eye 
is provided for each object in parenthesis. The scale bar in each 
image represents 2" . Figures ^] 1-^31 are available in the online 
version of the Journal. 



ternal M31 extinction. For each test field, we insert 46 
artificial clusters into individual raw (FLT) images us- 
ing new functionality developed for the DOLPHOT pho- 
tometry package. The clusters are randomly positioned 
within each image to prevent search bias (as would result 
from a regular grid pattern) , though we ensure that clus- 
ters are well-separated within the image so that they do 
not overlap or interfere with the photometry of other ar- 
tificial clusters. We drizzle the resulting images together 
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Fig. 5. — Same as Figure^] but for possible clusters listed in 
Table[2]with 2.0 < S hy _ eye < 2.5. Figures [5] l-[5] 12 are available 
in the online version of the Journal. 

to create searchable images in the optical F475W and 
F814W passbands. In all, we created three fields of arti- 
ficial clusters for each of the four full bricks in the Year 
1 dataset, resulting in a total of ~550 simulated objects. 
In addition to the completeness tests that follow, these 
artificial clusters are also used for quality assu ranc e of 
our photometric results, as discussed in Section |42) 

After creating the artificial cluster images, we identify 
clusters in the same way we searched the original im- 
ages, with all eight team members searching each field 
and rating the reliability of the clusters. By comparing 
the resulting cluster identifications to the full list of in- 
serted clusters, we can estimate the completeness limits 
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validation analysis in Section[0] Photometry zeropoints 
for the ACS and WFC3 cameras were obtained from the 
STScI webpagj^l and are listed in Table [3] All photom- 
etry is presented in the Vega magnitude system using 
the native HST passbands; we do not perform passband 
conversions. 



Fig. 6. — Cluster completeness derived from artificial cluster tests. Left: Completeness in the two optical filters, derived from 3 pc artificial 
clusters inserted in B09, B15, and B21. Center: Variation in completeness as a function of brick membership, derived from R e ff — 3 pc 
artificial clusters. Right: Variation in completeness as a function of R e ff, derived from artificial clusters inserted in B09, B15, and B21. 

of the sample based on various cluster input parameters; 
the results of this analysis are shown in Fig. [6j The left 
panel shows the best estimate of the sample's characteris- 
tic cluster recovery fraction, computed using a subsample 
of 110 "typical" simulated clusters with R e ff of 3 pc dis- 
tributed within the outer three bricks (B09, B15, & B21). 
We estimate a 50% completeness limit of m^AnhW ~21.2. 
Considering distance and Galactic foreground dus t red- 
dening (E(B -V) = 0.062; |Schlegel et aT||1998| ), this 
translates to an absolute magnitude limit 01 Mfa75w ~ 
-3.5. 

The true completeness of the sample as a whole, how- 
ever, is a complicated function of cluster luminosity, size, 
and location within M31. Using the full complement of 
simulated clusters, we explore the dependence of com- 
pleteness on the latter two cluster properties. The middle 
panel of Fig. [6] shows that variations in galactocentric po- 
sition result in a ~1 magnitude difference in completeness 
limits. Detection limits are brighter for clusters located 
in the bulge-dominated inner galaxy due to luminous, 
crowded background fields. Cluster size also plays a role, 
in which increasing size at constant luminosity results in 
reduced surface brightness and lower detection efficiency. 
The right panel of Fig. [6] shows that there is a ~1 magni- 
tude difference in completeness when comparing clusters 
with R e ff of 1 versus 7 pc. To account for these multiple 
dependencies, we will undertake a larger and more rig- 
orous set of completeness tests in future work to better 
characterize the subtleties of the completeness function. 



4. INTEGRATED PHOTOMETRY 

We use aperture photometry to measure the six band 
integrated fluxes of each cluster. Aperture photometry 
consists of two main analysis tasks: defining an aperture 
(center and size) and determining the local background 
flux lev el. We build upon photom etry procedures used in 
M3Tby [Barmby fc Huchra| fl200ip and |Krienke fc Hodge 
( 2007 ) , with several refinements related to the assessment 
ol local background levels and accounting for light that 
lies outside the photometric aperture. A detailed de- 
scription of photometr ic an alysis procedures and results 
is provided in Section |4.1[ followed by artificial cluster 



4.1. Aperture Photometry Procedure 

The first step in aperture photometry is to define an 
aperture center. Centers are estimated by centroiding on 
a F475W image that is smoothed using a 0.3" (6 pixel) 
FWHM Gaussian kernel. For certain clusters, particu- 
larly low luminosity objects whose light is dominated by 
a small number of bright sources, the flux-weighted po- 
sitions determined by the automated procedure do not 
always accurately reflect the cluster center. For this rea- 
son, we visually inspect central positions and manually 
adjust incorrect determinations. 

We adopt a photometric aperture size that provides 
the largest signal-to-noise ratio for flux measurement by 
enclosing a maximum amount of cluster light while in- 
cluding as little background light as possible. The clus- 
ters considered in this study vary by a factor of ~10 in 
radius, and consequently, the chosen photometric aper- 
ture radii vary by the same amount. We define circular 
apertures using growth curve analysis to determine an 
appropriate radius. The aperture limit is defined at the 
radius where the cluster profile drops below the level of 
the noise in the background, equivalent to the point at 
which the curve of growth turns over and the increase in 
cumulative flux as a function of radius stops. An illustra- 
tive example of a cluster image and growth curve is pro- 
vided in Fig. [7j The aperture radii (R ap ) are determined 
by visual inspection of the growth curves for each cluster 
and are reported with the photometric results. While it 
would be preferable to adopt an algorithmic approach for 
defining R ap , the relatively noisy character of the local 
background significantly complicates automated deter- 
minations. We perform aperture definition and growth 

19 http://www.stsci.edu/hst/acs/analysis/zeropoints and 
http :77www^sTscT™. edu/hst/wf c3/phot_zp_ffin] 
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curve analysis on the F475W image, which provides the 
best combination of signal-to-noise and contrast between 
cluster and field populations for a wide range of cluster 
ages. Apertures of the same angular size are used for 
the five other images. Photometric aperture radii range 
between 0.5" and 6" for the Year 1 sample. 

Aperture photometry depends greatly on the deter- 
mination of the background flux level. Following tra- 
ditional photometric terminology, we also refer to non- 
cluster background light as the "sky" or the "sky back- 
ground" . For the PHAT cluster sample, the sky is made 
up of two components: individual resolved stars and un- 
resolved light. Traditionally, the background is deter- 
mined using the mode of sky region pixel values. How- 
ever, the fact that resolved stars are a true component 
of the background light calls for an alternate statistical 
treatment. 

We define ten annular sky regions that encircle the 
photometric aperture, extending radially from 1.2 x R ap 
to 3.4 x R ap . Each of the annuli are equal in area to the 
photometric aperture, to accurately measure the disper- 
sion of the sky background based on an equal number of 
pixel samples. Figure [7] shows an example of the aperture 
layout, where we denote the inner and outer extent of the 
sky measurement annuli. Next, we calculate the total in- 
tegrated flux within each of the sky regions. We perform 
iterative 2a rejection on these ten fluxes, thereby exclud- 
ing regions that contain bright stars or other objects not 
representative of the typical sky background. We adopt 
the mean of the non-rejected sky fluxes as the sky back- 
ground value, and the standard deviation of these fluxes 
as a measure of the uncertainty in the sky background 
determination. We propagate the uncertainty from the 
sky level determination into the overall cluster photom- 
etry by adding it in quadrature with the cluster's flux 
measurement uncertainty. 

In agreement with the HKC cluster studies, the uncer- 
tainty in the sky background determination dominates 
the overall uncertainty in cluster magnitudes. When 
compared to previous ground-based cluster photometry 
in M31 or other HST-based extragalactic cluster photom- 
etry, our magnitude uncertainties appear larger. Our in- 
clusion of sky level uncertainty into the reported values 
account for these larger overall magnitude uncertainties. 
Uncertainty in the determination of cluster fluxes, com- 
parable to the errors reported in other cluster catalogs, 
never rise above 0.01 mag for any object in the Year 1 
catalog. 

The resulting magnitudes measured within R api here- 
after referred to as aperture magnitudes, are presented in 
Table [I] for the cluster sample and Table [2] for the possi- 
ble cluster sample. These aperture magnitudes represent 
high signal-to- noise, spatially matched measurements of 
cluster light, and are optimal for calculating cluster col- 
ors. We plot the photometric uncertainties as a function 
of magnitude for the cluster sample in Fig. [8j To as- 
sess relative reliability in the six passbands, we count 
the number of well-determined (a < 0.5 mag) photomet- 
ric measurements in each band, and record the results 
in Table [4j The F475W imaging provides the highest- 
quality measurements, as shown by that filter's small 
photometric uncertainties, followed by the F336W image 
which probes a similar wavelength regime. The F814W 
image has increased levels of photometric error due to the 



reduced contrast between cluster and field populations. 
The quality of the measurements is lower in the three 
remaining filters due to intrinsic wavelength-dependent 
limitations; the F275W measurements suffer from low 
signal-to-noise for all but the youngest, bluest clusters, 
while the F110W and F160W measurements suffer from 
high sky background levels and poor cluster-field con- 
trast. In addition to the problem of faint signal, we note 
that -10% of F275W magnitudes (and F336W magni- 
tudes, but at a lower level) are affected by cosmic ray 
artifacts. While a vast majority of these defects are ade- 
quately corrected for by our image processing, we caution 
that a few percent of the UV magnitudes might still be 
adversely affected. 

In addition to the aperture magnitudes, we also pro- 
vide estimates of the effective radius (R e ff, equivalentlv, 
the half-light radius) for each cluster in Tables [I] and [2] 
These estimates are obtained by measuring the cluster 
flux profile and interpolating this curve to find the radius 
that contains half of the light within R ap . We use the 
F475W cluster light profiles to make these assessments, 
given their good d ata quality. We discuss the result- 
ing sizes in Section |7.2| and use them to calculate aper- 
ture corrections in the next section. We recognize that 
the R e f f estimates will systematically underestimate the 
true cluster sizes due to the fact that they are calculated 
using aperture magnitudes that fall short of measuring 
the full luminosity of the cluster. However, we expect 
the impact of this underestimation to be small due to 
the steepness of the cluster light profile. 

4.1.1. Aperture Corrections 

The aperture magnitudes presented above measure a 
majority of the cluster light. However, these values do 
not account for light that lies below the noise level of the 
sky background, beyond the limits of our photometric 
aperture in the faint outer wings of the cluster profile. 
We correct for this missing flux by using R e ff determi- 
nations to approximate the cluster's luminosity profile 
shape, then extrapolate this profile past the limits of our 
photometric aperture to make an estimate of the cluster's 
total light. The magnitude difference calculated between 
the original aperture magnitudes and this total light es- 
timate are equivalent to an aperture correction. 

To make this profile-based extrapolation, we require a 
cluster profil e shape and a normalization for that profile. 
We adopt a King] ( [l962| ) profile, assuming a concentra- 
tion (c = Rudai I Rcore) that matches the characteristic 
profile of PHAT clusters, as obtained during preliminary 
cluster profile fitting (c = 7). Next, we use a cluster's 
R e ff to set the radial scaling of that characteristic King 
profile. Finally, we normalize the scaled profile such 
that the integrated flux within R ap matches the cluster's 
aperture magnitude measurement. Once normalized, we 
calculate the fraction of flux that lies outside R ap and 
transform this value into an aperture correction in mag- 
nitudes, which is presented in Tables [T] and [2] for all clus- 
ters. These corrections may be applied to the aperture 
magnitudes to obtain total magnitudes^] These total 
magnitudes are optimal for the estimation of total cluster 
magnitudes and luminosities. While the aperture correc- 

20 Total Magnitude = Aperture Magnitude + Aperture Correc- 
tion 
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Fig. 7. — Cluster aperture layout and growth curve for PC362. Left: F475W image of cluster with aperture locations overplotted. 
The green circle denotes the photometric aperture. The red circles denote the inner and outer edges of the annulus used for background 
determination; this region is subdivided into ten equal-area annular subregions to estimate the variance in the background. Right: The 
flux profile (top) and cumulative growth curve (bottom) for the cluster; vertical red and green lines correspond to the radii plotted in the 
left panel. 
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Fig. 8. — Photometric errors for integrated cluster photometry in each of the six PHAT passbands. The scale of the y-axis varies between 
panels. A standard reference at a = 0.5 mag is provided for comparison. 



tions were derived in the F475W passband, they may be 
used for all filters under the simplifying assumption of 
flat radial color profiles in the outer parts of the cluster. 

The amplitude of the aperture corrections are pre- 
sented in Fig. [9j Over the sample of clusters, the cor- 
rections vary from 0.0 to -0.6 mag, with a median cor- 
rection of -0.1 mag. These corrections are negligible for 
the brightest clusters, where a majority of the light is 
detectable above the noise level of the sky background. 
The corrections become larger for fainter clusters, due to 
their low cluster-to- field flux contrast. The simplifying 
assumption of a universal cluster profile shape that varies 
only as a function of R e f f provides suitable accuracy for 
this correction, as sh own by artificial cluster tests that 
follow in Section |3~2) Full cluster profile and structural 
parameter fitting is currently underway (Fouesneau et al. 
2012, in prep.), and these results could be used to further 



refine these profile extrapolations and improve upon our 
estimates of total cluster light in future work. However, 
we expect little overall improvement in photometry as a 
result of increased aperture correction precision because 
these corrections are comparable in size to the amplitude 
of the photometric uncertainties for most clusters. 

4.2. Artificial Cluster Photometry Experiments 

We use artificial clusters to assess the uncertainties and 
biases associated with our photometry procedures. For 
these tests, we only consider clusters that lie in the outer 
three bricks, have R e ff of 3 pc, and were dete cted as part 
of our artificial cluster search in Section |3.2| Analysis of 
this sample of 76 artificial clusters provides an evalua- 
tion of photometric accuracy for "typical" clusters in the 
Year 1 sample. We process the simulated objects using 
photometric procedures identical to those described in 
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Fig. 9. — Aperture corrections, used for converting aperture mag- 
nitudes to total magnitudes. The corrections are derived from clus- 
ter profile extrapolation. 
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Fig. 10. — F475W magnitude differences between input and re- 
covered magnitudes for a subset of simulated clusters inserted into 
B09, B15, and B21, adopting a R e ff of 3 pc. Top: Differences be- 
tween the measured total magnitudes and the input magnitudes. 
Bottom: Differences between the measured aperture magnitudes 
and the input magnitudes. The lowest outlier in both panels (cir- 
cled) differs due to a bright (F475W~19) nearby field star that 
overlaps the cluster. 



the previous section. 

The results of these tests are shown in Fig. [lOj The 
artificial cluster tests show good agreement (AF475W 
<0.1 mag) between the input and output magnitudes for 
bright clusters (F475W < 19). For less luminous clusters, 
aperture magnitudes are fainter than input magnitudes 
by up to ~0.4 mag at F475W ~ 21.5. However, total 
magnitudes are more successful in recovering input mag- 
nitudes accurately, showing a smaller faintward bias of 
~0.2 mag at F475W ~ 21.5. This bias represents a ~0.5cr 
deviation when compared to the ~0.4 mag scatter in the 



photometry. Finally, we acknowledge that these pho- 
tometric experiments are slightly idealized, for instance 
due to our consideration of only a single family of cluster 
profile shapes. However, this testing confirms that our 
photometry techniques provide accurate assessments of 
cluster properties and their associated uncertainties. 

5. COMPARISON TO EXISTING M31 CLUSTER STUDIES 

As discussed in the introduction, there is a long his- 
tory of stellar cluster studies in the Andromeda galaxy. 
Decades of effort have produced a wealth of knowledge 
on this topic. To place our findings in context, in this 
section we cross-reference our cluster identifications with 
existing catalogs, allowing us to reference previous work 
on the same objects and compare the results of our clus- 
ter analysis to existing catalogs. 

To compile a list of known clusters located within 
the Year 1 footprint, we bega n by searching the Re- 
vised Bologna Catalog (RBC; |Galleti et aL]|2004| last 



updated 2009 December to v4aJJ! This excellent re- 
source has aggregated all known cluster identifications 
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Huxor et al. 2008 


Caldwell et al.|2009p. Other than the KBC, we searched 



cent RBC revision (IVansevicius et al.|2009| |Peacock et al. 



2010} |Fan et al ||2010p . Our search of the HK(J pro 
duced addit ional o b j ect s for cross-matching, however no 
new objects were recovered from the three other catalogs. 
The Peacock et al. (2010) catalog is composed solely of 
edits and reclassifications from an earlier ver sion of the 
RBC (v3.5). Similarly, the objects studied in |Fan et al.| 
( 2010[ ) are derived directly from v4.0 of the RBC. Finally. 



we find no overlap between any part of the PHAT foot- 
print (existing or plann ed coverage) and the so uthwest 
region of M31 studied by Vansevicius et aT| ( |2009|) . In the 
discussion that follows, we consider cluster classifications 
with respect to those provided in the RBC v4.0. 

The RBC and HKC contain a total of 146 published 
clusters, 32 cluster candidates, and 84 other non-cluster 
classifications (foreground stars, background galaxies, 
and HI I regions) that lie within the Year 1 PHAT sur- 
vey footprint. We cross-match these previously known 
objects with all identifications made as part of the Year 
1 search (cluster, possible clusters, and unlikely objects). 
Further, we re-examined the PHAT data at positions of 
previously cataloged objects that were not matched to 
PHAT candidates to ensure the complete re-analysis of 
all existing catalog entries that lie within the Year 1 foot- 
print. In total, we classify 132 of the previously known 
objects as clusters, 12 as possible clusters, and reject the 
118 remaining objects. We note that all confirmed Year 
1 PHAT clusters were identified independently as part of 
our blind cluster search. Tables [l] and [2] provide cross- 
m atched names of clust ers using the naming convention 
of |Barmby et al. (2000), consisting of the identifier from 
the Revised Bologna Catalog, followed by the identifier of 
the next most significant cluster catalog, where such ex- 
ists. Additional details concerning the comparison and 
reclassification of existing cluster identifications is pro- 
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Fig. 11. — Magnitude histogram comparing the Year 1 PHAT 
clusters (black) to objects from existing catalogs that fall within 
the Year 1 survey footprint (red). The inset shows a breakdown 
of the previously known clusters into objects that were identified 
using ground-based imaging (green) and those that were identified 
using HST imaging (blue). 

vided in Appendix [B] There we provide catalog-specific 
commentary on the validity of the previously published 
classifications, as well as object-by-object classification 
revisions for the RBC and the HKC. 

A comparison between ground-based and space-based 
M31 cluster catalogs reveals the importance of high spa- 
tial resolution imaging in cluster identification work. In 
Fig. [TT] we present a histogram of the number of clus- 
ter as a function of apparent magnitude for the Year 
1 cluster sample of 601 clusters and the 132 previously 
known objects confirmed as clusters in the PHAT data. 
It is immediately apparent from the plot that HST imag- 
ing has allowed the PHAT cluster survey to identify 
hundreds of low luminosity clusters that could not be 
identified from ground-based data. As shown in the in- 
set of Fig. [TT] the sample completeness associated with 
previously known clusters discovered from ground-based 
data (green dashed histogram) drops precipitously at 
mFA7bW > 18 {MpA7bW > —6.5). This fall-off in com- 
pleteness reflects the difficulty in differentiating between 
single unresolved stars and compact clusters in low reso- 
lution images. In contrast, high spatial resolution imag- 
ing from HST enables the identification of clusters 2-3 
magnitudes fainter than previous ground-based surveys. 
Within this fainter luminosity range, Fig. 11 shows an 
order of magnitude increase in the number of objects 
identified in the PHAT cluster catalog when compared 
to previous HST-based cluster survey work (HKC; blue 
dotted histogram in inset). This improvement results 
from the order of magnitude increase in spatial cover- 
age provided by PHAT when compared to the limited 
number of previous, targeted HST observations that fall 
within the Year 1 survey footprint (35 arcmin 2 versus 
390 arcmin 2 in PHAT Year 1). 



In addition to the catalog comparison presented here, 
in Appendix [C] we compare photometry results presented 
in this work to those of existing catalogs. This analysis 
acts as quality assurance for the photometry presented 
here, and provides the reader with an assessment of the 
inherent differences between the sets of photometric re- 
sults. 

6. YEAR 1 CLUSTERS: PHOTOMETRIC PROPERTIES 

The Year 1 cluster sample, derived from ~l/4 of the 
total expected PHAT survey data, provides the first 
glimpse of what can be expected from the full balance of 
the PHAT stellar cluster survey. In Section [5] we showed 
that our catalog represents a considerable increase in the 
number and diversity of clusters known in M31. Our ex- 
cellent HST-based imaging should lead to factor of >4 in- 
creases in the number of known clusters within the PHAT 
survey footprint. 

To obtain a better sense for the type of clusters we 
have identified in the Year 1 sample, we present color 
and magnitude distributions from the catalog photome- 
try. For this analysis, we select a subset of objects with 
well-determined photometry, where the uncertainties in 
the F336W, F475W, and F814W magnitudes are each 
less than 0.5 mag. This quality cut results in a sub- 
sample of 482 well-characterized objects that we use to 
explore the properties of the catalog. 

We plot a color-color diagram for clusters with well- 
determined photometry in Fig. 12 This diagram aids 



our ability to assess the cluster age distribution. To pro- 
vide reference points to guide the eye, we overlay the stel- 
lar evolution model predictions from the Padova group 
(Girardi et al. 2010| ). These model tracks assume solar 
metallicity and are reddened to account for foreground 
Galactic extinction (E(B - V) = 0.062; |Schlegel et all 
1998). On this plot, ages increase as we follow the evolu- 



tionary track from the upper left to the bottom right. In 
addition to the initial cluster selection based on photo- 
metric uncertainties, we define a second subset of clusters 
based on cluster luminosity. The red points denote the 
well-determined subsample's 92 most luminous clusters, 
with F475W < 19.5. 

The first conclusion we draw from Fig. [12] is that the 
Year 1 cluster sample includes a wide range of ages. Clus- 
ters populate the full length of the model evolutionary 
track, with a large number of objects populating an in- 
termediate age range (300 Myr to 3 Gyr). However, the 
mapping from position in the color-color diagram to age 
suffers from well-known degeneracies with extinction and 
metallicity. For example, old (12-14 Gyr) metal-poor 
([Fe/H] < —1.0) globular clusters inhabit the same posi- 
tion on this diagram as ~100 Myr old, solar metallicity 
clusters with Ay of ~1.5 mag. The green region in the 
right panel of the figure denotes the shared color-color 
region where this particular age-metallicity degeneracy 
exists. 

Second, we observe that the cluster color distribution 
is affected by the effects of stochastic stellar mass func- 
tion sampling in low-mass (< 10 4 Mq ) clusters. The 
red region in the right panel of Fig. [12] highlights color 
outliers that have anomalous red F475W-F814W colors. 
While young clusters (<50 Myr) suffering from large 
amounts of dust attenuation (Ay >2 mag) can popu- 
late this region of the diagram, these red colors are more 
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Fig. 12. — Left: Color-color diagram for the 482 Year 1 clusters with well-constrained (a < 0.5 mag) photometry in the F336W, F475W, 
and F814W passbands. A subset of the 92 most luminous clusters (F475W < 19.5) are highlighted in red. Padova SSP models for Solar 
metallicity, reddened to account for Galactic foreground extinction, are plotted as a cyan line for reference. Cyan points act as age indicators, 
spaced at 0.5 dex increments beginning at 10 Myr. Median error bars for the bright (red) and faint (black) cluster samples are displayed 
in the lower left corner. Right: A smoothed, grayscale version of the color-color diagram, with SSP models shown again as reference. The 
red box denotes the location of anomalous red supergiant clusters (RSGCs), which contain luminous evolved stars that strongly affect the 
cluster's integrated colors. The green box denotes the parameter space populated by metal poor globular clusters (MPGCs). The blue line 
denotes the modeled location of the cluster main sequence (CMS), reddened by Ay =0.4 mag to match the color distribution of the Year 1 
cluster sample. The CMS represents the color-color sequence populated by low mass clusters that host no evolved stars due to stochastic 
sampling of the cluster's stellar mass function. 

frequently caused by the presence of a small number of 
bright evolved supergiant stars that can bias integrated 
cluster colors. We d iscuss this behavior in greater de- 
tail in Section |7.3.1| In addition to the red outliers, 
stochastic effects can also cause integrated cluster colors 
to appear bluer than model predictions. The fluctuation 
in the small number of evolved stars sometimes results 
in the complete absence of supergiant cluster members, 
meaning that the cluster's integrated light is emitted ex- 
clusively by main sequence stars. Such clusters fall onto 
a linear sequence in color-color space we refer to as the 
cluster main sequence. We highlight this feature in blue 
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in the right panel of Fig. 12 
Next, we plot a cluster color- magnitude diagram in 



Fig. [13] to assess the cluster mass range probed by the 
Year 1 sample. As in Fig. [l2j we plot foreground red- 
dened, solar metallicity Padova stellar models for refer- 
ence. The PHAT clusters span ~8 mag in F475W lumi- 
nosity, translating to >3 orders of magnitude in cluster 
mass. This range indicates that the PHAT cluster sam- 
ple hosts a wide variety of clusters, spanning systems 
that contain hundreds of solar masses up to those with a 
million solar masses. The most luminous clusters, how- 
ever, all appear to have red colors (F336W-F814W ~ 
2.5), forming a vertical sequence of objects on the right 
side of the plot. This results from the fact that most 
massive (> 10 5 M ) clusters in the Year 1 sample are 
old globular clusters. 

Figures [12] and 13 show that the sample of clusters as- 
sembled from the PHAT dataset provide a top-to-bottom 
assessment of the M31 cluster population. Few datasets 
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Fig. 13. — Color-magnitude plot for the 482 Year 1 clusters with 
well-constrained (<r < 0.5 mag) photometry in F336W, F475W, 
and F814W passbands. Padova SSP models for three cluster 
masses (10 5 , 10 4 , 10 3 Mq) at Solar metallicity are plotted for 
reference, with age indicators spaced at 0.5 dex increments begin- 
ning at 10 Myr. Characteristic median error bars are shown on the 
right side of the plot, where each point describes the uncertainties 
of points within a 2 mag bin in luminosity. 
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have the ability to sample objects across a variety of 
stages in cluster evolution over such a large, uninter- 
rupted mass range. This diversity makes the PHAT clus- 
ter sample a valuable tool to better understand the for- 
mation and dissolution of star clusters. 

7. DISCUSSION 

7.1. Luminosity Functions 

Luminosity functions are a basic, model independent 
measure of a stellar cluster population. Up to this point, 
observations of spiral galaxies have provided largely con- 
sistent results, where luminosity functions are well fit 
by a pow er- law (N oc L~ aL ) with indices of about —2 o r 
steeper ( [Larsen|200^|Gieles|2Q10{IChandar et al.|2010b| ). 
In addition, there is growing evidence of steeper slopes 
at the brightest cluster luminosities (My < —9; |Whit- 



th 



Imore et al.|fl999 Gieles et al.||2006| iHaas et al.|[2008), 



suggesting the possibility of a Schechter-like cluster mass 
funct ion with an e xponential truncation at the high-mass 
end flLarsen||2QQ9l ). 

These previous studies have generally focused on the 
bright end of the luminosity function, fitting clusters with 
luminosities greater than ~ 2 x 10 4 L (equivalent to 
My < —6). In contrast, the faint end of luminosity func- 
tion has relatively poor constraints due to detection and 
catalog completeness difficulties encountered by previous 
extragalactic cluster studies. However, due to our abil- 
ity to identify low luminosity clusters in PHAT survey 
data, we can use the Year 1 cluster sample to probe the 
shape of the luminosity function down to limits only pre- 
viously accessible in the Magellanic Clouds . These new 
constraints are interesting because of the sensitivity of 
the faint end luminosity function to cluster disruption 
flLarsen||2QQ9l ). 

To assess the shape of the PHAT Year 1 luminosity 
functions in each of the six filter passbands, we consider 
two different samples of objects for fitting: all 601 clus- 
ters, and a subsample of 534 disk clusters that excludes 
objects from the bulge-dominated Brick 1. The cluster 
population within Brick 1 is dominated by old, massive 
globular clusters associated with the galaxy's bulge com- 
ponent, in addition to the significantly brighter complete- 
ness limits in this region. The disk sample represents a 
set of predominately younger objects that are likely more 
homogeneous in terms of formation and destruction his- 
tory, simplifying the interpretation of its resulting lumi- 
nosity functions. We narrow the samples further to ex- 
clude clusters with uncertain photometry (a > 0.5 mag), 
performing the quality cut individually for each pass- 
band. The fraction of objects in each filter that meet 
this quality requirement is listed in Table [4j Finally, 
we convert aperture corrected, total magnitudes for the 
selected clusters to absolute magnitudes, correcting for 
the M31 distance modulus and Galactic foreground red- 
dening. We note that no correction for dust attenuation 
within M31 has been applied; we will derive these correc- 
tion factors on an object-by-object basis as part of future 
age and mass fitting analysis (Fouesneau et al. 2012, in 
prep.). 

To characterize the power law slope of the luminos- 
ity function, we perform a linear fit to constrain the 
slope of log dN/dM (the logarithm of the number of 
clusters per magnitude) as a function of absolute mag- 



nitude. For each passband, we bin the clusters using 
variable-sized magnitude bins to ensure fair weighting 
of the data, fol lowing the suggestion of D'Agosti no fc 
Step hens] fl!986|) (and more rece ntly in |Mafz Apellaniz &; 
Ubeda||2005| |Haas et al.||2008| to group the data such 
that each bin represents an equal number of clusters, N. 
We choose N = 15, but find the results are insensitive 
to the particular number chosen. We fit to datapoints 
brighter than absolute magnitude completeness limits: 
-4.0 for F275W, F336W, and F475W; -5.0 for F814W; 
—6.0 for F110W and F160W. These limits are conserva- 
tive; they are equivalent to > 80% completeness for all 
bricks, as found in Section 3.2 We convert the resulting 
magnitude-based slopes to equivalent luminosity func- 
tion slopes (oll) an d report these as our primary results. 
The luminosity functions and their associated power-law 
fits for the complete and disk-only cluster samples are 
presented in Figs. [14] and [15] the values of the fitted 
slopes are listed in Table [5] We separate the results into 
two groups to isolate the F275W, F110W, and F160W 
passbands, which carry the potential for larger system- 
atic uncertainties due to the smaller fraction (~50-75%) 
of clusters with well-determined photometry available at 
these wavelengths. As such, we will focus our subse- 
quent discussion on results from the F336W, F475W, 
and F814W passbands in Fig. [14] 

Considering the complete sample of Year 1 clusters, we 
find luminosity function power law slopes in the F336W, 
F475W, and F814W passbands that are all flatter than 
—2. While these measurements agree with the general 
trend of flatter slopes me asured at fainter l uminosities 
(e.g., in NGC45 and M51; |Mora et al.||20Q7j 

)es result ri 



Haas et al. 
Tom the in- 



V v,. ft ., i,, unu L , | Mora 

2008), we suggest these flat slopes 
elusion of a relatively large number of luminous, evolved 
globular clusters associated with Brick 1 and included in 
the complete cluster sample. We find that although clus- 
ters in Brick 1 make up ~15% of the complete sample by 
number, ~50% of sample members with rriF475W < 19 
(Mf475\v i$ —5.5) are Brick 1 clusters with spectroscopi - 
cally derived ages of >10 Gyr ( [Caldwell et al.|2009[|2011[ ). 
Globular clusters are known to follow a Gaussian-like lu- 
minosity function shape with a peak magnitude around 
My of —7 and —8, depending on factors such as metal - 



licity and age (e.g., Harris 1991 Barmby et al. 2001) 



The dissimilar luminosity function shape of the globular 
cluster subpopulation significantly influences the slope of 
the overall complete Year 1 luminosity function. 

To determine the luminosity function shape that is 
characteristic of the dominant, young cluster popula- 
tion found in M31, we examine the fitting results for 
the disk-only sample of clusters, presented in the bot- 
tom row of Fig. [14] We find uniformly steeper slopes 
(a L < -1.9) for the F336W, F475W, and F814W disk- 
only luminosity functions when compared to the com- 
plete Year 1 sample results. In addition to the closer 
agreement between these slopes and the canonical —2 
power law, we also observe a significant trend where the 
luminosity functions at bluer wavelengths are flatter than 
at redder wavelengths. This behavior can be explained 
by wa velength and age-dependent cluster mass-to-light 
ratios flGielespMo] ). 

To visualize the sensitivity of the faint end luminos- 
ity function slope to differences in cluster dissolution, we 
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Fig. 14. — Luminosity functions for the F336W (left), F475W (center), and F814W (right) passbands. Top panels display the luminosity 
functions for clusters drawn from the complete Year 1 sample, and the bottom plots show the luminosity functions considering only 
disk clusters that lie outside the bulge (i.e., excluding B01 clusters). Only clusters with well-constrained (a < 0.5 mag) photometry are 
included in the analysis. Clusters have been corrected for Galactic foreground dust attenuation, but not attenuation within M31. Data is 
plotted using variable binning, such that each point represent an equal number of clusters (N=15). Vertical error bars represent Poisson 
uncertainties due to the number of objects per bin, and horizontal error bars represent the magnitude range of each bin. The luminosity 
function slope is fit down to completeness limits denoted by vertical dotted lines at -4.0 for F336W and F475W, and -5.0 for F814W. 



compare our F475W luminosity functions to shapes pre- 



dicted by four canonical dissolution models in Fig. [16 
The analytica l luminosity functions were calculated by 



Larsen ( 2009 ) under the assumption of a constant cluster 
formation history, an underlying Schechter mass func- 
tion (with cutoff of 2xl0 5 M©), and continuous, non- 
stochastic sampling of the stellar initial mass fun ctional 
In the case of mass-dependent cluster dissolution, Larsen 
(2009) finds that the slope of the luminosity function 
should flatten to > —2 at My > —8 or —9, while sce- 
narios with mass-independent dissolution (and the case 
of zero dissolution) result in slopes that always remain 
steeper than —2. 

The complete Year 1 cluster sample appears to be 
better described by a mass-dependent disruption model, 
while the steeper disk cluster luminosity function appears 
to agree with the shape of the mass-independent or the no 
disruption case. However, due to the degenerate effects of 
underlying cluster age and mass distribution differences 

22 We note that in terms of luminosity function modeling, we 
find little evidence for any noticeable differences when stochastic 
sampling of a cluster's stellar initial mass function is incorporated, 
validating the present model comparison. 



and variations in cluster disruption, such simple interpre- 
tation is not possible. We opt to defer further analysis 
and interpretation until we obtain robust age, mass, and 
extinction determinations for individual clusters. This 
crucial information will allow for the separation of the 
degenerate effects that influence the shape of the lumi- 
nosity function and better constrain the characteristics 
of the underlying cluster population. 

7.2. Cluster Sizes 

Cluster size is a fundamental parameter that provides 
constraints on the dynamical state of stellar systems. 
While many interesting correlations involving cluster 
sizes require the determination of other cluster param- 
eters (e.g., mass-radius and age-radius relations), from 
sizes alone we can compare the overall properties of the 
cluster sample to those in other galaxies and look for in- 
teresting sub-populations of clusters within the sample. 

We plot R e ff as a function of F475W magnitude in 
Fig. [l7| In this parameter space, the sample separates 
well into two groups: Brick 1 clusters that are predom- 
inately old (~10 Gyr), massive globular clusters as soci- 
ated with the galaxy bulge as discussed in Section |7.1 
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Fig. 15. Same as Figure (14] but for the F275W (left), F110W (center), and F160W (right) passbands. Completeness limits are -4.0 for 
F275W and -6.0 for F110W and F160W. 



and clusters from the remaining bricks that are associ- 
ated with Andromeda's star forming disk. Overall, the 
Brick 1 clusters are more compact and luminous than the 
disk clusters, and show a trend in which R e ff decreases 
with decreasing F475W luminosity. Along with the clus- 
ter data, we plot a line of constant surface brightness an- 
chored at F475W of 21.2 and a R e ff of 3 pc, representing 
the 50% catalog completeness limit as determined in Sec- 
tion [372] The distribution of clusters with respect to this 
line suggests that, to first order, catalog completeness 
roughly follows a constant surface brightness cutoff. 

The disk cluster sample includes objects with a broad 
range of ages and suffers from significant completeness 
truncations, making its magnitude-radius distribution 
difficult to interpret. On the other hand, the Brick 1 clus- 
ters are similarly aged and are generally much brighter 
than the completeness limit. This allows us to conclude 
that the moderate trend in which fainter Brick 1 clus- 
ters are more compact is real. We do not currently know 
the physical origin of this trend, but plan to fol l ow up 
using ages and metallicities fro m Ca ldwell et al. (2011) 
and masses from Strader et al. ( 2011| ). 

The left panel ol Fig. |1 (^presents the principal re- 
sult of our size analysis, which shows that the overall 
size distribution of the Year 1 PHAT cluster sample can 
be described approximately as a log- normal distribution, 
where the median R e ff value is 1.68 pc (0.44"). This me- 



dian R e ff is similar to, but smaller than va lues found in 
other galaxies. In M83, [Bastian et al.|(|2011| ) find a over- 
all median val ue of ^2.4 pc lor a sample ol mostly young 



clusters, while iMcLaughlin & van der Marel 



( [2005] ) find a 

median value ol ^3.3 pc lor a sample ol Milky Wa y glob- 
ular clusters. Additionally, iBarmby et al. (2007) find a 



median value of ~2.5 pc for a sample ol previously known 
M31 globular clusters. 

We can explain the differences in median cluster sizes 
by considering how the composition of the Year 1 PHAT 
cluster sample differs from the M83, Milky Way, and ex- 
isting M31 samples. One possible reason for the me- 
dian mismatch with the Milky Way and M31 globular 
cluster samples is that the Year 1 PHAT sample is a 
heterogeneous mix of mainly young clusters, while the 
globular samples are uniformly much older. However, 
an age-based explanation is ruled out because the old, 
globular-dominated Brick 1 subsample has a similar, but 
still smaller median R e ff of 1.69 pc. The R e ff distri- 
bution for the Brick 1 subsample is shown in the right 



panel of Fig. 18 



An explanation for the median differences becomes 
clear when we consider the galactocentric radii of the 
existing globular cluster samples. A correlation between 
R e ff and galactocentric radius has been show to exist 
in ma ny galaxies such that inner clusters are more com- 
pact flBarmby et al.||2007[). When we limit the Milky 
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Fig. 16. 

shows the luminosity function constructed using all 
a subset of Year 1 disk clusters that lie outside the bulge-dominated B01. The outcom es of four different cluster dissolution models are 
compared: no disruption, two models of mass-dependent disruption as parameterized by |Lamers et al^| (|2005) for two dissolution timescales 
(£4; characteristic destruction timescale for a 10 4 Mq cluster), and mass-independent disruption in which 90% of clusters are destroyed 
within every logarithmic age interval. Models are normalized to the data at M.jp4^5w — — 5. Although the samples appear to prefer 
particular dissolution models, the effects that complex underlying cluster age and mass distributions impart on the luminosity functions 
do not allow us to draw any conclusions about cluster disruption at this time. 

Way and previous M31 cluster samples to include only 
objects with galactocentric radii less than 3 kpc, which 
roughly matches the radial extend of the Brick 1 sample, 
the median R e ff for these two samples drops to ~2.0 
pc. These lower medians are in better agreement with 
the Brick 1 subsample, easing the tension between the 
values derived for the two globular cluster samples and 
for the PHAT clusters. 

To explain the median R e ff difference between the 
young M83 cluster sample and the Year 1 PHAT sam- 
ple, instead of age and galactic position, we consider how 
completeness limits influence the shape of the R e ff dis- 
tribution. In terms of numbers, the Year 1 cluster sample 
is dominated by compact, low luminosity objects that lie 
just above our detection threshold. At these faint mag- 
nitudes, surface brightness limits prevent us from detect- 
ing clusters across the full range of possible R e ff values. 
As a result, the overall cluster sample is biased towards 
compact objects, driving the median towards smaller val- 
ues. However, if we only consider clusters with F475W 
< 20.5, the biasing effects of surface brightness limits are 
greatly reduced. The R e ff distribution for this luminous 
subsample is shown in the right panel of Fig. [THJ With- 
out the faint, compact portion of the sample, the median 
cluster size rises to 2.32 pc, which agrees well with the 
observed value obtained for clusters in M83. In conclu- 
sion, cluster sizes for the Year 1 sample are consistent 
with results derived by previous works in similar spiral 
galaxies once sample completeness and composition are 
taken into account. 
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Fig. 17. — Reff as a function of F475W magnitude. The old, 
globular cluster dominated population in B01 (red solid circles) 
occupy a different region of parameter space than the younger disk 
clusters that dominate the remaining outer bricks (blue open cir- 
cles). The black dashed line denotes an approximate 50% com- 
pleteness limit, corresponding to a uniform surface brightness limit 
consistent with a F475W cutoff of 21.2 at an R e ff of 3 pc. 
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Fig. 18. — Histograms of R e ff derived from F475W data. Left: Distribution of R e ff fc> r the full Year 1 cluster sample. Median R e ff is 
1.68 pc (0.44"). Right: Distribution of R e ff for the B01 (red) and bright disk (green) subsamples. Median R e ff for the two subsamples 
are 1.69 pc (0.44") and 2.32 pc (0.61"), respectively. 



7.3. Objects of Interest 
7.3.1. Red Supergiant Clusters 

We identified a subsample of objects with unusually 
red F475W-F814W colors in Fig. [12} Within the selec- 
tion box, we identify ~15 clusters that show extreme 
integrated colors. Upon inspection, we find that many 
of these clusters appear to host luminous red supergiant 
(RSG) stars, which explain their anomalous integrated 
colors. Two examples of this class of cluster, PC57 and 



PC1127, are shown in Fig. 19 



Red supergiant stars are massive (8 to 25 M ) stars 
that have left the main sequence and are traversing the 
top of the CMD during the helium burning phase of stel- 
lar evolution. Individual luminous RSG stars have the 
ability to bias the integrated colors of low mass (< 10 4 
Mq) clusters because stochastic sampling of the stellar 
mass function in low mass systems can cause the small 
number of evolved stars to vary by factors of a few (low 
mass clusters typically host between and 4 evolved 
stars). This variation creates wild fluctuations in the to- 
tal luminosity and integrated color of the cluster, which 
depend on the particular number of evolved stars pro- 
gressing through this relatively short evolutionary phase. 

The existence of these color outliers confirms that stel- 
lar population modeling of low mass clusters must ac- 
count for stochastic sampling of the stellar initial mass 
function to obtain accurate age and m ass determinations 
(see e.g., jFouesneau fc Lancon] 2010). Models that as- 
sume a fully sampled stellar initial mass function cannot 
reproduce objects with integrated colors as red as those 
shown by RSG clusters. Our future cluster analysis will 
benefit from the use of cluster models that account for 
stochastic fluctuations in integrated light caused by the 
RSGs. 

More importantly, our ability to resolve individual 
stars in clusters provides us a number of important ben- 
efits with respect to the RSG clusters. Resolved star 



photometry enables us to use CMD fitting techniques to 
obtain cluster parameter determinations, allowing us to 
avoid the use of biased integrated measurements all to- 
gether. In addition, once we obtain age, mass, and atten- 
uation information for the clusters, we can use this infor- 
mation to tag individual RSGs and use these constraints 
to improve calibration of late stage stellar evolution for 
massive stars at high metallicity. 

7.3.2. Massive Clusters 

Globular clusters, with old characteristic ages (~10 
Gyr) and large masses (> 10 5 Mq), have long been a 
target of study in M31. Although the Year 1 cluster 
sample is numerically dominated by low mass clusters, it 
contains many massive clusters as well. While we iden- 
tified few new globular clusters as part of the Year 1 
search, our high spatial resolution imaging has enabled 
us to confirm or eliminate a large number of globular 
cluster candidates from existing catalogs. In addition to 
the 63 out of 71 highly-ranked clusters we confirm from 
the RBC, we affirm 12 and reject 19 possible candidates. 
The PHAT survey also enables detailed analysis of these 
objects by means of resolved star photometry for cluster 
members. We plan to place better constraints on red gi- 
ant branch (RGB) and horizontal branch (HB) morphol- 
ogy as a function of metallicity through these detailed 
studies. 

In addition to the old globular clusters, we have also 
newly identified a number of intermediate mass (~ 10 4 
M ), intermediate age (~1 to 3 Gyr) clusters. These 
objects are interesting targets f or study be cause of the 
relative rarity in the Milky Way ( |Friel|1995| ), where most 
similarly aged Galactic clusters are less massive. While 
a small number of these intermediate mass, interm ediate 
age clusters are alrea dy known to exist in M31 ( |Cald-| 
well et al.|2QQ9 , 2011), increasing the sample size for this 
class ol object should help to better understand their 
origin and evolution. Investigation of these objects will 
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Fig. 19. — Resolved star color-magnitude diagrams and color optical (F475W+F814W) image cutouts for two candidate red supergiant 
clusters. These clusters have unusually red F475W-F814W colors due to the presence of luminous evolved stars that strongly bias integrated 
light measurements. 



comple ment the study of younge r (<1 Gyr) massive clus - 
ters by|Fusi Pecci et al.| fl2005| ), |Caldwell etH] ( |2009| ), 
and JPerma et al.| ( |2010| T 

The most remarkable massive cluster in the Year 1 
catalog is the cluster SK142C (PC1156). This object 
is the sample's most discrepant color outlier, appearing 
at the redward extreme of all cluster photometric mea- 
surements, with a F475W-F814W color of -3.5. This 
cluster appears to be an intermediate mass, intermedi- 
ate age, highly-reddened cluster. Along with cutout im- 
ages showing the cluster's heavily reddened appearance, 
we plot the cluster's resolved star CMD in Fig. 20 The 
CMD shows an extended, highly-reddened RGB, indica- 
tive of an intermediate cluster age. We perform a vi sual 
comparison to so lar metallicity Padova isochrones ( |Gi-| 



rardi et al. 2010) and find that the CMD is consistent 
with an age ol —1-3 Gyr and Ay of 2.25-2.75 mag, as- 
suming large uncertainties associated with visual "chi- 
by-eye" isochrone fitting. Using this combination of age 
and reddening, we use the cluster's F160W luminosity 
to estimate a cluster mass of — 1.5 x 10 5 M . Due to 
its massive nature and its somewhat intermediate age, 
this cluster is a rare and interesting ob ject akin to newly 



discovered objects in the Milky Way ( Fig er et al.||2006 
Davies et aL|2011a|b ). 



8. SUMMARY & FUTURE WORK 

In this paper, we presented the first installment of 
the PHAT stellar cluster catalog. We introduced the 
Year 1 cluster sample, consisting of 601 clusters iden- 
tified through a visual search of high spatial resolution 
HST imaging. The PHAT cluster sample significantly 
increases the number of clusters known in M31; the cat- 
alog presented here represents more than a factor of four 
increase in the number of known clusters within the cur- 
rent survey area. We presented a basic assessment of the 
cluster sample, including positional, size, and photomet- 
ric information. 

We have shown that the PHAT cluster sample hosts a 



large range in ages and masses. This wide-ranging sam- 
ple provides the opportunity for a top-to-bottom study of 
stellar cluster evolution processes. The PHAT sample in- 
cludes low-luminosity objects missed in studies of distant 
external galaxies, while covering an uninterrupted range 
of cluster masses, unlike Milky Way cluster samples. 

When the survey is complete, the PHAT cluster cat- 
alog will be among the largest and most comprehensive 
surveys of star clusters in any galaxy. This work presents 
results derived from the first 25% of the survey data; we 
estimate that the final sample will include ~2500 clus- 
ters. Over the duration of the survey, we plan to peri- 
odically publish updates to the catalog to include new 
clusters and to revise object classifications as we gather 
additional information about the sample. Age and mass 
determination analysis for the Year 1 sample is currently 
underway (Fouesneau et al. 2012, in prep.), which will 
provide the means to explore a host of different topics, 
including the cluster mass function and cluster dissolu- 
tion. The Year 1 cluster sample will also be the basis for 
analysis of structural parameters, resolved star content, 
and integrated spectroscopy, but we also look forward to 
future studies enabled by the complete, four year PHAT 
cluster catalog. 
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Fig. 20. — Left: Image cutouts showing highly-reddened, intermediate age cluster SK142C (PC1156). Right: The resolved star color- 
magnitude diagram for SK142C shows a reddened RGB, indicative of an intermediate age population. We overplot a Padova stellar 
isochrone with a log(Age) of 9.35 (~2 Gyr), Ay of 2.5 mag, and solar metallicity to show that the cluster is consistent with these parameter 
values. 



Facilities: HST (ACS, WFC3). 



APPENDIX 

BACKGROUND GALAXY CATALOG 

As introduced in Section |3.1[ we present a catalog of 
objects visually identified as background galaxies in Ta- 
ble [6] Along with positional information, we present a 
rough assessment of the galaxy size. 

COMMENTARY ON EXISTING CLUSTER CATALOGS 

We performed a comparison of object classifications 
from existing cluster catalogs to those in the PHAT Year 
1 cluster catalog in Section [5] Here we provide comments 
on individual catalog comparisons, including detailed 
classification statistics and object-by-object classification 
revisions. We discuss ind ividual classification statistics 
for the RBC (IGalleti etaf 



2004 



it jKim et all (|2007| ) , [Cald- 
well et al.| ( |2fl(J9J |2Ulip 7[Feacock e t al.| (|2Ul( J|h and three 



ol the lour HKC studies ( |Krienke & Ho dge 2007; Hodge 
2010) no objects lrom |Krienke & Hodge| 2008 
e Year 1 PHAT footprint). Classification re- 



et al._ 
lie within t 



include 71 clusters, 36 cluster candidates, and 80 non- 
cluster objects. After reanalyzing these classifications 
using the by-eye ranking results and following up on 
other objects that were not selected as PHAT candi- 
dates, our revised object classifications show excellent 
agreement with the original RBC results. Individual ob- 
ject revisions are provided in Table [8j but the bulk of the 
changes are due to our ability to determine the nature 
of the RBC's candidate objects (RBC Flag = 2). Of the 
36 candidates, we confirm 12 objects as clusters in the 
PHAT catalog, reject 23 candidates (including two orig- 
inally controversial cases, one duplicate entry, and one 
object not found in PHAT imaging), and transfer the fi- 
nal candidate to the PHAT possible cluster sample. In 
terms of catalog accuracy, only five "confirmed" clusters 
(RBC Flag = 1) were confide ntly revised to no n-cluster 
objects, four of which were |Kim et al.| ( |2007[ ) objects 
(see additional discussion below). Overall, the high level 
of consistency between the two catalogs gives us great 
confidence in galaxy- wide accuracy of the RBC cluster 
s ample 



Kim et al. 



Kim et al. 



suits broken down by object type are presented in Table 

m 

The Revised Bologna Catalog — Our search for RBC 
objects that lie within the Year 1 PHAT imaging foot- 
print produces a catalog of 187 objects. These objects 



(2001) — The catalog content of 
(2007]) is also accounted for in the results of the RBC 
cross-comparison, but the resulting cluster confirmation 
and rejection statistics from t his work is worth indiv idual 



mention. As p ointed out in [Caldwell et al. (2009) and 
|Peacock et"aL] fl2010p , the |Kimet al | ( |2UU7| ) catalog has 
a high level ol contamination, namely from misclassified 
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stellar sour ces. Consid ering all three cluster quality cate- 
gories from Kim et al. (2007) together, we find 38 objects 
from this catalog that fall within the Year 1 PHAT foot- 
print. Of these, only 37% of those objects appear in the 
PHAT catalog as confirmed or possible clusters. Even 
when taking the quality rankings into consideration, we 
find that (50%, 27%, 38%) of the (A, B, C) objects ap- 
p ear in the catalogs from our present work. 
\Caldwell et al\ {200^ \2011\ ) — These two papers re- 



vised classifications from v3.4 of the RBC (many of 
which were adopted in v4.0) using low-resolution spec- 
tra, ground-based imaging from the Local Group Galaxy 
Survey (LGGS; [Massey et al.|2006[ ), and HST-based im- 
ages. As this catalog has been incorporated into the lat- 
est version of the RBC, discussion of individual object 
classifications are included in Table [8] When comp aring 
our current PHAT classification with those from I C aid- 1 
we find good agreement between 



well et al. 



(2009, 



>]|2oTT]), 

We revif 



the two works. We revise classifications for 6 clusters out 
o f 183 that fall w i thin t he Year 1 footprint. 

Peacock et al \201Cty — The clusters considered by 
Peacock et al. "(^0TI?)are included in the RBC compari- 
son discussed above, and classifications for each of these 
objects is already provided in Table [8} Many of the re- 
classifications performed in|Peacock et al. (2010[), which 



was based on v3.5 of the RBC, have proven to be accu- 



rate (inclu ding the stellar classification of numerous Kim 
|et al.|200"7| c andidates). Overall, w e find good consistency 
between the |Peacock et al.| ( |2010[ ) catalog and our PHAT 
classifications. 

Hodge- Krienke Catalogs — Similar data and methodol- 
ogy (including common authors in the case of P.H.) leads 
to strong consistency between the HKC and the work 
presented here. We find 54 clusters from the HKC that 
match entries in the PHAT cluster catalog and 8 addi- 
tional objects that match possible clusters, out of a total 
of 75 objects that lie within the boundaries of the PHAT 
Year 1 footprint. Unlike in the case of ground-based cat- 
alogs, the 13 HKC objects not recovered as part of the 
PHAT cluster search mostly fall into the category of as- 
terisms. In other words, these objects were interpreted 
by the PHAT search team as uncertain or unlikely cluster 
candidates, and omitted from the final catalog presented 
in this work. This subset of objects shows that for the 
least luminous, least massive clusters, defining the differ- 
ence between an object that is potentially a bound stel- 
lar system and simply a chance collection of stars seen in 
projection becomes a difficult and subjective task. Table 
[9] presents proposed revisions to the HKC. In closing, we 
note that we do not compare t o clus ter identifications 
presented in Williams fc Hodge] ( |2001[ ) because classifi- 
cations in this work were superseded by the HKC. 

COMPARISON TO EXISTING CLUSTER PHOTOMETRY 

We provide a comparison between PHAT cluste r pho 
tometry and t he results of [Peacock et al. J 2010|) 
et al l fl2010D , |Hodge et al] ( |2009| ), ano^ 
( |2010 ) . We compare photometry lor these selected stud- 
ies because they provide magnitudes derived in similar 
photometric passbands as part of a uniform analysis. 
We choose not to compare to photometry provided by 
the Revised Bologna Catalog due to the heterogenous 
nature of their unified photometric measurements (see 



Fan 



lodge et al. 



Galleti et al. 2004 , for details) . We not e that both Pea- 
^2010D and |Fan et al.| ( |2010|) find go od agree- 



cock et al. 



ment betwe en photometric results and that of |Barmbyl 
et al. (2000), upon which the photometry of the Revised 
Bologna Catalog is based. 



The optical aperture photometry of Pea cock et al. 



(2010) was derived from SDSS ugriz imaging. We 



compare photometric results for 77 clusters found in 
both datasets. We calculate transformations to convert 
from AB-based g and z-band photometry to Vega-based 
F475 W and F814W mag nitudes using transformations 
from |Girardi et aL] ( |2008| ), assuming a single-age, 1 Gyr 
solar metallicity population that represents a median 
value for the age-dependent transformation: 

F475W = g AB + (0.06 ± 0.02), (CI) 

F8UW = iab - (0.54 ± 0.10). (C2) 

The resulting photometric comparison is presented in 
the top row of Fig. |2l] We find that the PHAT mag- 
nitudes are fainter than those of Peacock et al. (2010) 
by ~0.25 mag. This offset can be explained by th e dif- 
ference in apertur e sizes, as seen in Fig. 22 The Pea- 
cock et al. (2010) apertures have radii that are twice 



as large on average, and up to four times as large as 
those used in this work. In the case of the brightest 
comparison cluster, B127-G185 (PC1425), the size of the 
aperture (10") was large enough to include the nearby 
cluster NB89 (PC1426), explaining the -0.5 mag differ- 
ence for this object. In a comparison of F475W-F814W 
versus (g — i) optical colors, the agreement is quite good 
for brighter clusters (F475 W <17) with in c reasin g scat- 
ter and a redward bias of Peacock et al.| (|2010[ ) colors 
for fainter clusters. These differences ar e likely caused 
by fie ld contamination within the large | Peacock et al.| 
(2010) apertures, and by blending in their low-resolution 
(1-2" seeing) ground-based imaging. In summary, we 
find systematic differences be tween the photometr y de- 
rived in this work and that oflPeacock et al.l (|2010|), but 



the overall agreement is adequate and remaining differ- 
ences can be readily explained by the effects of image 
resolution and aperture size differences. 

In another rece nt compilati o n of g round-based aper- 
ture photometry, Fan et al. ( 2010| ) derived Johnson- 
Kron-Cousins UBVR 1 magnitudes from LGGS imaging 
flMassey et aL||2006| >. The Vega-based PHAT F336W 
and F814W magnitude s are nearly equi valent to the U 
and I magnitudes from |Fan et al.| (2010| ), allowing us to 
directly compare photometric results lor 68 clusters com- 
mon to both samples. We present the results in the sec- 
ond row of Fig. [2lJ s howing good agr eement between the 
two datasets. While [Fan et al.| ( [2010| ) magnitudes appear 
slightly brighter for faint clusters, trie overall agr e ement 
is improved with respect to the 



Peacock et al. 



(2010) 



results. One explanation for the improvement between 
ground-based and HST-based results could trace back to 
the superior image quality of the LGGS data when com- 
pared to the SDSS imaging of |Peacock et al. (2010). The 
average seeing for the LGGS data is while varying 
between 1 -2" for the SPS S imaging. This results in the 
ability for Fan et al. (2010) to reduce their aperture sizes, 
such that they are only a factor of ~ 1.5 bigger than those 
used by PHAT. 



The photometry of Hodge et al. ( 2009 ) and Hodge et al. 
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Fig. 21. — Comparison of P HAT photom etric results with those of |Peacock et al.| ( |2010| (top row), [Fan et al.| ( |2010| ) (second row),[fiodge 
[eTaL] |2009t (third row), and Hod ge et al.] |2010| (bottom row). 

(I2010D allows us to compare two sets of HST-based pho- 
Out of the four M31 cluster catalog 



tometric results, 
papers that compose the HKC, we choose two where the 
tabulated photometry is provided in native, Vega-based 
HST filter magnitudes. Although the data used in these 
studies are similar to those used in this work, there are 
fundamental differences in the photometry techniques 
employed. As opposed to our growth curve method of 
aperture size determination, the HKC adopt an isopho- 
tal aperture determination, such that apertures extend 
out to a chosen surface brightness limit. Further, the 
methods of sky background determination differ. Our 
sky levels are determined using local sky background es- 
timates measured in annular rings around the photomet- 
ric aperture, while the HKC uses samplings of the sky 
background taken across the HST image. Due to these 
differences in methodology and the lack of any correction 



to account for total cluster light, we expect the compar- 
ison to show some scatter and a systematic offset such 
that the PHAT photometry i s brighter. 

For the |Hodge et al. ( |2QQ9[ ) comparison, we use 18 ob- 
jects common to both aatasets and present the results 
in the third row of Fig. 21 We compare our F475W 



photometry to magnitudes measured in a similar F450W 
passband, using a transformation derived in the same 



way as Eqs. |C1| and |C2 



F475W = F450W - (0.11 ± 0.10). 



(C3) 



There are several significant (A >0.5 mag) outliers in the 
F475W magnitude comparison, in which the photometry 
from this work is brighter by up to ~1 mag. To check the 
validity of our measurements, we performed photometry 
for the three most discrepant out liers on the WFPC2 
images used by |Hodge et al.| ( |2009| ) , employing photom- 
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etry procedures used in this work. We found resulting 
magnitudes that agree within ~0.1 mag of our original 
PHAT photometry values. This result indicates that the 
magnitude offsets are likely due to differences in photo- 
metric technique, specifically related to sky background 
determination. 
The bottom row of Fig. |2l] presents the comparison of 



Hodge et al. (2010) photometry for 13 common objects. 
We find better o verall agreement w ith this dataset when 



compared to the|Hodge et al. (2009) results. We compare 



F814W magnitudes (no filter transformations required) 
and F336W-F814W colors, and find a scatter of -0.2 
mag between the two datasets. The comparison of pho- 
to metry between our current PHAT study and the work 



of|Hodge et al. (2009, 2010) shows good overall consis- 



tency. Systematic biases and a small number of signifi- 
cant outliers likely stem from differences in measurement 
technique. 



R Qp (PHAT) [arcsec] 



Fig. 22. — C omparison of photometric aperture radii used in 
|Peacock et al.| ( |2010| and in this work. The solid line represents 
a 1-to-l relation, the dotted line a 2-to-l relation, and the dashed 
line a 4-to-l relation. 
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TABLE 1 
PHAT Year 1 Cluster Catalog 



PC ID 

^by — eye 


RA (J2000) 
riap { ) 


Dec (J2000) 

K eff { ) 


F275W 

rpoi A\\T 


a 

(7 


F336W 

r 1 1U VV 


a 

(7 


F475W 

r lOUVV 


a 

(7 


ApCor a 
Alternate Name 


1 

1.00 


11.638274 
1.57 


42.193887 
0.68 


17.54 
18.34 


0.02 
0.04 


17.68 
17.95 


0.03 
0.07 


18.82 
17.73 


0.02 
0.12 


-0.12 
HodgelO-85 


2 

1.43 


11.637139 
2.51 


42.209936 
0.64 


15.60 
17.23 


0.02 
0.04 


15.91 
17.39 


0.02 
0.19 


17.33 
17.30 


0.01 
0.35 


-0.00 
HodgelO-84 


20 
1.00 


11.630550 
1.10 


42.200631 
0.48 


19.83 


0.07 


22.55 
19.16 


0.25 
0.23 


21.97 
18.63 


0.63 
0.39 


-0.13 


21 
1.57 


11.631591 
1.00 


42.199991 
0.52 


19.60 
20.98 


0.02 
0.28 


19.82 
22.88 


0.01 
1.80 


20.89 


0.10 


-0.21 


22 
1.71 


11.630849 
0.75 


42.201656 
0.46 


23.38 
21.39 


0.36 
0.67 


22.52 
22.47 


0.11 
1.72 


22.22 


0.06 


-0.31 


Note. — 
regarding it 


Tablets published 
s form and content. 


in its entirety in the electro: 
The sample presented here c 


nic edition of the Astrophysical Journal. A p 
onsists of objects classified as clusters, wi 


th S by _ 


shown here for guidance 
eye < 2 -0. 


a Aperture 


Corrections are proA 


rided such that r, 


a Total — m 


Aperture 


+ ApCor. 











TABLE 2 

PHAT Year 1 Possible Cluster Catalog 



PC ID 


RA (J2000) 


Dec (J2000) 


F275W 


a 


F336W 


a 


F475W 


a 


ApCor a 


^by — eye 


Rap (") 


R eff (") 


F814W 


a 


FllOW 


a 


F160W 


a 


Alternate Name 


4 


11.664001 


42.192391 


19.37 


0.06 


19.43 


0.05 


20.38 


0.02 


-0.09 


2.43 


2.10 


0.83 


19.80 


0.14 


19.42 


0.22 


18.63 


0.10 




6 


11.668817 


42.197452 


21.50 


1.01 


20.93 


0.13 


21.86 


0.29 


-0.01 


2.29 


1.13 


0.32 


21.97 


1.72 


21.57 


0.69 


21.71 


3.33 




7 


11.675393 


42.193555 






25.59 


2.19 


22.63 


0.14 


-0.07 


2.29 


0.80 


0.29 


20.96 


0.13 


20.38 


0.11 


19.73 


0.08 




11 


11.648244 


42.227958 


25.42 


1.35 


23.98 


0.24 


22.56 


0.15 


-0.08 


2.43 


1.00 


0.38 


20.12 


0.12 


18.77 


0.03 


17.82 


0.07 




12 


11.661360 


42.190571 


21.39 


0.16 


21.02 


0.05 


21.38 


0.04 


-0.02 


2.00 


1.05 


0.31 


20.74 


0.06 


20.82 


0.13 


21.39 


1.90 





Note. — Table|2jis published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. The sample presented here consists of objects classified as possible clusters, with 2.0 < S^y _ eve < 2.5. 



Aperture Corrections are provided such that m-Total = 771 Aperture ApCor. 



TABLE 3 
Photometric Zeropoints 



Passband 


Vegamag Zeropoint 


F275W 


22.65 


F336W 


23.46 


F475W 


26.16 


F814W 


25.52 


FllOW 


26.07 


F160W 


24.70 



TABLE 4 

Passband Photometric Quality Comparison for Cluster Sample 



Passband N (Valid Measurements) N (Well-determined Measurements) 



F275W 


552 (91.8%) 


447 (74.4%) 


F336W 


590 (98.2%) 


566 (94.2%) 


F475W 


600 (99.8%) 


597 (99.3%) 


F814W 


593 (98.7%) 


514 (85.5%) 


FllOW 


518 (86.2%) 


358 (59.6%) 


F160W 


472 (78.5%) 


313 (52.1%) 


Note. — Valid me 


asurements denote magnitudes that r< 


ssult from positive fluxes (signal 


measured above sky 1< 


3vel) and suffer no other failures (e.g., 


image artifacts). Well-determined 


measurements denote 


magnitudes where er < 0.5 mag. 
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TABLE 5 
Luminosity Function Fits 



Passband 


cm I, (All Clusters) 


a l (Disk Clusters) 


F275W 


-1.72 ±0.07 


-1.76 ±0.08 


F336W 


-1.73 ±0.06 


-1.90 ±0.08 


F475W 


-1.79 ±0.06 


-2.06 ±0.09 


F814W 


-1.74 ±0.05 


-2.18 ±0.09 


F110W 


-1.68 ±0.06 


-2.20 ±0.11 


F160W 


-1.56 ±0.05 


-1.82 ±0.08 



TABLE 6 

Background Galaxy Catalog 



ID 


Brick 


RA 


Dec 


R 






(J2000) 


(J2000) 


(arcsec) 


1 


21 


11.677181 


42.192120 


1.41 


2 


21 


11.644169 


42.197941 


1.20 


3 


21 


11.647937 


42.202160 


1.37 


4 


21 


11.685157 


42.218735 


0.97 


5 


21 


11.616282 


42.196312 


2.05 



Note. — Table \6l is published in its entirety in the elec- 
tronic edition of the Astro-physical Journal. A portion is shown 
here for guidance regarding its form and content. 



TABLE 7 

Summary of Existing Cluster Catalog Classifications and Revisions 



Catalog Clusters Candidates Galaxies HII Regions Stars Other a Total 



RBC 


78 


(71) 


4 


(34) 


3 


(3) 


1 


(3) 


98 


(76) 


3 


(0) 


187 


Kim (Total) 


12 


(10) 


2 


(28) 


2 


(0) 





(0) 


22 


(0) 





(0) 


38 


Kim A 


4 


(10) 


1 


(0) 


2 


(0) 





(0) 


3 


(0) 





(0) 


10 


Kim B 


4 


(0) 





(15) 





(0) 





(0) 


11 


(0) 





(0) 


15 


Kim C 


4 


(0) 


1 


(13) 





(0) 





(0) 


8 


(0) 





(0) 


13 


Caldwell 


68 


(73) 


1 


(0) 





(0) 


1 


(0) 


113 


(110) 





(0) 


183 


Peacock 


77 


(58) 


4 


(25) 


2 


(2) 





(2) 


41 


(40) 


3 


(0) 


127 


HKC 


54 


(75) 


8 


(0) 





(0) 





(0) 





(0) 


13 


(0) 


75 



Note. — In each column, the first number represents the number of objects per category after 
PHAT reclassification and the second number (in parenthesizes) repres ents the number of original 
classifications from the published catalog. The A classification from |Kim et al.| (2007) maps to a 
cluster classification in the RBC, while B and C map to candidate classifications. 
a The "other" classification signifies clust ers that were not reco vered by the PHAT search or were 
duplicate catalog entries in the RBC and Pea cock et al.| ( [2010| ) catalogs. In the case of the HKC, 
the "other" classification signifies objects that were deemed non-cluster asterisms. 



TABLE 8 

Revised Bologna Catalog Revisions 



Cluster Name 


PC ID 


New Flag 


Old Flag 


Comments 


BH16 


1381 


1 


2 




B523 


1383 


1 


2 




SK118C 


641 


1 


2 




SK134C 


1349 


1 


2 




M028 


544 


1 


2 





Note. — Table [8] is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. 
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TABLE 9 
Hodge-Krienke Catalog Revisions 



Cluster Name 


PC ID 


PHAT Classification 


WH13 


1119 


Possible Cluster 


WH18 


1089 


Possible Cluster 


KHM31-195 


1282 


Possible Cluster 


KHM31-241 


964 


Possible Cluster 


Hodge09-57 


977 


Possible Cluster 



Note. — Table |9] is published in its entirety in 
the electronic edition of the Astrophysical Jour- 
nal. A portion is shown here for guidance regard- 
ing its form and content. 



